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Abstract. We slightly extend the fluctuation theorem obtained in [21] for sums 
of generators, considering continuous-time Markov chains on a finite state space 
whose underlying graph has multiple edges and no loop. This extended frame 
is suited when analyzing chemical systems. As simple corollary we derive by a 
different method the fluctuation theorem of D. Andrieux and P. Gaspard for the 
fluxes along the chords associated to a fundamental set of oriented cycles [3] . 

We associate to each random trajectory an oriented cycle on the graph and 
we decompose it in terms of a basis of oriented cycles. We prove a fluctuation 
theorem for the coefficients in this decomposition. The resulting fluctuation the- 
orem involves the cycle affinities, which in many real systems correspond to the 
macroscopic forces. In addition, the above decomposition is useful when analyz- 
ing the large deviations of additive functionals of the Markov chain. As example 
of application, in a very general context we derive a fluctuation relation for the 
mechanical and chemical currents of a molecular motor moving along a periodic 
filament. 

Keywords: nonequilibrium steady state, thermodynamic force, affinity, large de- 
viation, generating function, oriented cycle, spanning tree, fluctuation theorem, 
molecular motor. 

1. Introduction 

Out-of-equilibrium systems are common in daily life. Examples are mechanical 
systems in contact with thermal reservoirs as well reacting systems in contact with 
particle reservoirs generating particle fluxes through differences of chemical poten- 
tials. Considering chaotic dynamical systems of statistical mechanics, Gallavotti 
and Cohen |1U| have discovered a symmetry relation (now taking their names) for 
the large deviation functional of the average entropy creation rate. This relation is 
also called fluctuation theorem. The same result has then been derived for Markov 
stochastic processes by Kurchan |19| . Maes [25] and Lebowitz and Spohn [24] . An 
extensive rigorous treatment including further developments is given in the book of 
Jiang et al. [15] (see also references therein together with [6]). Near to equilibrium, 
the fluctuation theorem generalized to systems with several currents yields the On- 
sager's symmetry and the usual Green-Kubo's formulas for transport coefficients 
|11[ [24"] . In this sense, it can be thought as their generalization far from equilibrium. 
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The initial investigation of the fluctuation theorem has referred to models coming 
from statistical mechanics, and in the last years a proper analysis has been developed 
for chemical and biochemical systems by Andrieux and Gaspard (see for example 
[121 E El O HI [5] and references therein) . The analysis for a model of molecular motor 
along a periodic filament with two chemical states has also been done by Lacoste, 
Lau, Mallick [2(71 [2T| 1221 123] . 

The approach of Andrieux and Gaspard is much inspired by the network theory 
of out-of-equilibrium systems. This theory has been developed by Hill [13] and 
Schnakenberg [30] with a special attention to biochemical systems. Mathematically 
one is interested to the stochastic dynamics of a continuous time random walk on a 
finite connected graph Q with multiple edges and no loop, where the jumps along an 
edge have positive probability rate in both directions. For such a random walk the 
stationary distribution /x exists and is unique. To each oriented bond £ in Q from 
the state a to the state a', one associates a microscopic affinity defined as 



where A; (01,(72) denotes the probability rate for a jump from a\ to 02- Note that 
detailed balance in reversible systems simply corresponds to the fact that all micro- 
scopic affinities are zero. The affinity associated to a given oriented cycle C in Q is 
then defined as the sum of the microscopic affinities of the oriented bonds forming 
the cycle: A{C) = X^eC It is simple to check that the affinity of a cycle remains 

the same if one replaces the above microscopic affinity A(£) by ln(fc(0, a')/k(a',a)). 

The network theory presented in [30] is based on the fact that to each unoriented 
spanning (maximal) tree on the graph Q one can associate in a canonical way a 
family of oriented cycles, called fundamental set of oriented cycles, being a sort of 
basis of the space of all oriented cycles. Each oriented cycle in a fundamental set 
contains only one oriented edge (called chord) that does not belong to the spanning 
tree when disregarding orientation. In [30] Schnakenberg has observed that, for 
several models of real systems, the macroscopic (mechanical or thermodynamical) 
forces keeping the system out-of-equilibrium are the affinities of the oriented cycles 
of some fundamental set in Q (note that more oriented cycles can have the same 
affinity). In [3] this situation is called Schnakenberg condition. We point out (see 
[301 13]) that usually thermodynamic forces are not encoded in a single jump rate. 

In P|, starting from a spanning tree, Andrieux and Gaspard have derived a fluc- 
tuation theorem for the currents along the chords, whose conjugate variables are 
the affinities associated to the oriented cycles in the fundamental set. Under the 
Schnakenberg condition, this fluctuation theorem can be restated in terms of the 
macroscopic forces keeping the system out-of-equilibrium instead of cycle affinities. 

Let us explain our theoretical contribution in this direction. First, we recover 
the above results of Andrieux and Gaspard by a different approach. Indeed, the 
fluctuation theorem presented for sums of several generators in Section 2.3 of [24] 
is not suited for random walks on graphs with multiple edges, but as we show its 
simple proof can be easily adapted to the present case. This slightly extended result 
immediately leads to the fluctuation theorem for the family of currents along all 




li{cr)k{(j, a') 



(1.1) 
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oriented edges, and therefore to the fluctuation theorem for the currents referred 
to a fundamental set of oriented cycles. We point out that the fundamental sets 
of oriented cycles, canonically associated to the spanning trees of the graph Q, do 
not cover all possible bases of the cycle space. We then extend the network theory 
presented in [30] and prove a fluctuation theorem for the "currents" referred to a 
generic basis of oriented cycles. When dealing with a generic basis, chords disappear 
and the definition of the variables conjugated to the affinities of the oriented cycles 
in the basis is more algebraic. We call these conjugate variables generalized currents. 
Their definition needs to associate to each random trajectory an oriented cycle, which 
must then be decomposed in the basis. The generalized currents are the coefficients 
in this decomposition. Their definition is essentially algebraic and differs from the 
derived chain discussed in [15] [Chapters 1,2]. In the derived chain method, for each 
oriented cycle C one counts the number of times Wf(C) that the cycle has been 
realized by the trajectory up to a given time t (there is nothing algebraic): in [15] 
the authors show that this quantity, rescaled by t, converges a.s. to some circulation 
number lo{C) and show that the entropy production rate in the steady state can 
be expressed in terms of the circulation numbers. For a discussion on fluctuation 
theorems satisfied by uJt(C) we refer to [5]. 

The above cycle decomposition of the trajectory is particularly useful when look- 
ing for a fluctuation theorem of additive functionals of the Markov chain. As ex- 
ample of application in this direction, we consider a very general discrete model for 
a molecular motor moving along a periodic filament under the action of an exter- 
nal force /, transforming the chemical energy from ATP hydrolysis to mechanical 
work. The forces keeping the system out-of-equilibrium are the external force / 
and the chemical potential difference Afi associated to ATP hydrolysis. In the last 
years, much attention have been devoted to the thermodynamics of small systems 
supported by very fast technological improvements (see for example [U [29] and ref- 
erences therein). Molecular motors are proteins working as motors at the nanoscale, 
with very interesting thermodynamic aspects, which have been much investigated 
both from a theoretical and an experimental viewpoint |14[ [IB] \T7\ [28] . The fluc- 
tuation theorem for molecular motors with two or three chemical states have been 
proved in [T] and [201 \2l\ 122] I23j. The small number of chemical states allows a de- 
tailed analysis based on matrix computations, leading to more information than the 
Gallavotti-Cohen-type symmetry for large deviations. In the general case, matrix 
computations become not reasonable. On the other hand, the developed theory for 
fluctuation theorems associated to cycle decompositions allows to easily prove in full 
generality the fluctuation relation obtained in |21[ [22] for only two chemical states. 
More precisely, calling xt the position of the molecular motor along the filament and 
zt the ATP consumption, it holds 



where, roughly speaking, the function #(A,7) is characterized by the identity 



A, 7 e R, 



(1.2) 



t > 1 
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and where f(3 and x t denote the dimensionless quantities ff3 x [lm] and ^/[lm] 
([lm] being the length unit). Of course, j3 = 1/kT where k is the Boltzmann's 
constant and T the absolute temperature. 

1.1. Outline of the paper. In Section [2] we extend the fluctuation theorem for 
sums of generators (see Theorem I2.1j) . In Section [3] we discuss some fluctuation 
theorems (Fact [3T3l and Theorem 13. 5ft related to cycle decompositions (fundamental 
set of oriented cycles and basis, respectively). All proofs are postponed to Sections 
El [3 Finally, in Section |4] we discuss the biochemical application of Theorem 13.51 
leading to relation (11. 2\ . 

2. Fluctuation Theorem for sums of Markov generators 

As starting point, we extend the fluctuation theorem obtained in [24J [Section 2.3]. 
As in [24] we consider a continuous-time Markov chain (Xt)t>o on a finite state space 
S having Markov generator 

Cf(a)= Ka,a')(f(a')-f(a)), a € S . (2.1) 

We assume that the Markov chain is irreducible and that the transition rates satisfy 
the positivity relation 

k(cr,a') > <^ k(a',a) > 0. (2.2) 

The last assumption is very reasonable in physical systems, since for a real transition 
between physical states one expects that, at least with very low probability, the 
opposite transition is possible. We remark that, due to irreducibility and since <S is 
finite, there exists a unique stationary distribution of the Markov chain. 

We decompose the generator £ as a sum of m generators C^ 1 '. More precisely, we 
fix a family of numbers k^(a, a'), parameterized by the index i = 1, . . . , m and by 
the pairs (a, a') of distinct elements in S. We assume that: 

(Al) fcW(cr, a') > for all i: 1 < % < m and a ^ a' in S; 
(A2) YhLi k(i) ( a i a ') = o 1 ) for all a ^= a' in S; 

(A3) k^(a,a') > k^(a',a) > 0, for all i : 1 < i < m and a ^ a' in S. 

Trivially, the above family {A;W(cr, cr')|. , corresponds to a representation of C as 

sum of Markov generators, since we can write C = ^^Li £ where denotes the 
Markov generator on S defined as C^f(a) = ^2 a i :a '^ CT k ( - l \a,a')(f(a') — /(cr)). 

We point out that our assumptions (Al), (A2), (A3) coincide with the ones stated 
in [24] [Section 2.3], with the only exception that here we have dropped the additional 
assumption of [24] that, given a ^ a', there exists at most one index i such that 
k^(a,a') > 0. 

The above construction suggests to think of the Markov chain as a random walk 
on a connected graph Q with multiple edges and no loop: the vertexes of Q are 
given by the states of S, while two distinct states a, a' are linked by as many edges 
as the indexes i for which fcW (a,a'),k^ (a' , a) > 0, each edge is labeled by the 
corresponding index i. Then, k^ l \a,(j') represents the probability rate for a jump 
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along the i-labeled edge from a to a'. Such an interpretation of the Markov chain 
as random walk on a graph with multiple edges becomes particularly relevant in 
biochemical applications and is essential when considering thermodynamical forces 
and affinities [331 EJH Q21 E] • 

We formulate the fluctuation theorem in terms of the random walk on the graph 
Q presented above. By a trajectory up to time t, we mean the path (X s : s G [0, £]) 
together with the knowledge of the edges along which the walker has moved. Given 
i : 1 < i < m and states a ^ a', we define the weight 



w 



otherwise . 



To each trajectory up to time t, visiting the states o"o, (7%, . . . , a n (listed in chrono- 
logical order, n being a random integer) and jumping along the (unoriented) edges 
indexed respectively by iq, i±, . . . ,i n -i, we associate the functions W^\t), i : 1 < 
i < m, defined as 

n-l 

W®(t) :=Y, w{i) ^cj k+1 )5 ik=i . (2.4) 

fc=0 

Shortly, every time the walker jumps along an i-edge, the function W^> increases 
of the weight w^(a,a'), a and a' being the initial state and the final state of the 
jump. 

In what follows, given a distribution v on S we denote by ¥ u and E y the probability 
measure and the expectation w.r.t. the Markov chain Xt with knowledge of the 
crossed edges, starting with initial distribution v. If v = S a , we simply write P CT and 
E CT instead of P v and E w . 

We point out that for any Ai,...,A m G R, the moment generating function 



is finite. Indeed, all the weights w^(a,a') are finite, while 

by a simple coupling argument one gets that the random number of jumps in the 
time interval [0, t] is stochastically dominated by a suitable Poisson variable, which 
has finite moment generating function. 

Similarly to [23] [Section 2.3] (with the exception that here we consider arbitrarily 
initial distributions and we have dropped a technical assumption), the following 
holds: 

Theorem 2.1. Given a distribution v on S, for each Ai, . . . , A m G R the limit 

e(A 1 ,...,A m ) := lim -^lnE, f e -E£i^W] (2.5) 

t— >oc t L 

exists, is finite and does not depend on v. Moreover, the following fluctuation rela- 
tion is valid: 

e(Ai, . . . , A m ) = e(l - Ai, . . . , 1 - A m ) . (2.6) 
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The proof is given in Subsection 15.11 Due to the theory of large deviations (cf. 
[32j . Section 1.5 in [15] and references therein), the Legendre-Fenchel transform 

m 

I{z\,... ,z m ) := sup | e(Ai, . . . , \ m ) - ^ XjZj \ 

Ai,...,A m eK L ~[ ) 

is convex, lower-semicontinuous and non-negative. Moreover, it has compact level 
sets and satisfies inf zG ]Rm I(z) = 0. In addition, the random vector W^\t)/t satisfies 
w.r.t. P„ a large deviation principle with rate function / (independent from u). As 
discussed in [23], the fluctuation relation (|2.6p can be restated in terms of the rate 
function / as 

m 

I(z 1 ,...,z m )-I(-z 1 ,...,-z m ) = -J2zi, VzeM m . (2.7) 

i=l 

Physical implications of the above relations in the steady state are discussed in [24] . 

We point out that the above Theorem 12.11 leads to an infinite family of fluctuation 
relations ()2.6|) . parameterized by the decompositions C = YlTLi^ ■ Only some 
of them refer to relevant physical quantities. Some physically relevant cases are 
discussed in [23], the fundamental one corresponds to m = 1, in this case is 
related with the entropy production of the system. As explained in Section El by a 
suitable choice of the decomposition C = £ one immediately recovers from 

the above theorem the fluctuation theorems for currents obtained in [3]. 

3. Fluctuation theorem for currents 

The dynamical evolution of several physical and chemical systems is already de- 
scribed by a random walk on a finite graph, with multiple edges and no loop. Here, 
we can apply Theorem 12.11 to this context, obtaining an equivalent fact which is 
simply formulated in a different language. Moreover, we re-derive the fluctuation 
theorem for currents with respect to a fundamental set of oriented cycles obtained in 
[3] and consider a similar problem when working with more general sets of oriented 
cycles (Theorem 13.51 below) . 

We start now with a finite connected unoriented graph Q whose vertexes are given 
by the states in S. Q can have multiple edges but no loop, i.e. between two distinct 
vertexes a, a' there can be several edges, while there is no edge from a state to itself. 
We denote by £ Q (o = oriented) the set of edges of Q with orientation (each edge 
in Q has two possible orientations and therefore corresponds to two oriented edges 
in £ Q ). Given £ € £ , we denote by £ the edge obtained from £ by inverting its 
orientation. Moreover, we write £{ and if for the states in S such that t goes from 
£i to if (initial and final states). It is convenient to assign a canonical orientation to 
each unoriented edge in Q. We denote by £ c (c = canonical) the set of canonically 
oriented edges of Q. Note that each oriented edge in Q is given by i or i, for some 
£e£ c . 

We fix a family of positive numbers {k(£)}i e £ and consider the continuous time 
random walk Xt on Q, jumping along the edge £ (following the associated orientation) 
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with probability rate k(£). In particular, the Markov generator of Xt is given by 

Cf(a)= k W Wf) ~ /(*)] • (3-1) 



Since Q is a connected graph and since the constants k(£) are positive, X t belongs 
to the class of Markov chains on S introduced in Section [2j 
To each oriented edge i E £ , we associate the weight 

w(£) = \nk(£)/k(£) (3.2) 

and we define N#(t) as the number of times the random walk jumps along £ minus the 
number of times the random walk jumps along £, up to time t. In words, N^{t) is the 
flux along the oriented edge £. Considering the sum decomposition C = J2ie£ ^ 
where is the Markov generator associated only to the jumps along the oriented 
edges £ and £, it is trivial to check that Theorem 12.11 implies the following fact: 

Fact 3.1. For any initial distribution v and for any family {Xi : £ E £ c } of real 
numbers, the limit 

q({X e : £ E £ c }) := lim — lnE„ [ e -^ e£c >^(*)l (3.3) 
t— >oo t L 

exists, is finite and does not depend on u. Moreover, it holds 

q{{X e : £ E £ c }) = q{{w(£) - X e : £ E £ c }) . (3.4) 

We point out that, when v is the stationary distribution, the above result equals 
formula (84) in [3]. 

3.1. Fluctuation theorem for currents w.r.t. a fundamental set of oriented 
cycles. We recall here some basic concepts concerning the oriented cycles in a finite 
graph [131E01E]. [30] represents a very concise reference. 

An oriented cycle C'mQ can be described by a string of oriented edges (b\, . . . , 
bi E £ , such that the vertex in which bi enters equals the vertex from which bi+i 
exits (with the convention k + 1 = 1). We convey that the strings (pi, . . . , bk) and 
(bi, bi-i) identify the same oriented cycle C. When disregarding 

the orientation of C, we call it simply cycle. 

We fix a maximal tree (also called spanning tree) T on Q [30] , i.e. T is a 
unoriented subgraph of Q, containing all vertexes of Q and containing no cycle. The 
edges £ S £ c which do not belong to T when disregarding orientation are called 
chords of the maximal tree T. We enumerate the edges in £ c as l\, £2, ■ ■ ■ , £ n where 
£\,...,£ s are the chords associated to T. In particular, s is the number of chords of 
Q. It is simple to prove that s equals e — v + 1, where e is the number of edges of Q 
and v the number of vertexes. In particular, s does not depend on the choice of the 
maximal tree and will be referred in what follows as chord number. 

When adding to T a chord £i (disregarding its orientation), one obtains a graph 
containing a unique cycle. We give to this cycle the orientation induced by the 
chord £i and call C, the resulting oriented cycle. The set {C\, . . . ,C S } is called a 
fundamental set of oriented cycles. We now explain the origin of the name. Given 
an oriented cycle C and an oriented edge £ E £ Q , we denote by Si(C) the number of 
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times £ appears in C minus the number of times the reversed edge £ appears in C. 
Then 

s 

C = ^^.(C)C i . (3.5) 

3=1 

The meaning of the above identity is clarified by the following definition: 

Definition 3.2. Given oriented cycles C\, . . . , C^, C[, . . . ,C' r and given real numbers 
a% , . . . , Gtfc , a'i , . . . , a' r we set 

k r k r 

aid = ^2 a 'fii ^2 a i S ^ C i) = ^ a 'j S ^ C j) ^t££o- (3.6) 

i=l j=l i=l 3=1 

We point out that, given a fundamental set {C\, . . . ,C S }, the oriented edges ap- 
pearing in Cj are all distinct, that i and I cannot both appear in Cj and that, for 
each chord £i, it holds S^{Cj) = 5ij for 1 < i,j < s. 

Recall the definition of the weight ()3.2|) . Given a generic oriented cycle C we 
define its affinity as A(C) = Yli w {^)-> where the sum is among the edges t visited 
by the cycle with the proper orientation. Alternatively, one can set 

A(C) = St{C)w{l) . (3.7) 

We recall that in the above definition one can replace w(£) by the ratio between 
the local flux along £ and the local flux along £ w.r.t. the steady state or any 
other probability measure v on the vertex set, giving positive measure to each state, 
i.e. w{£) can be replaced by ]n\y(£i)k(£)/v(£f)k(£)]. This leads to an equivalent 
definition of cycle affinity. 

Recall the definition of the random variable Ng(t) given just after (|3.2|) . Ne(t) 
represents the flux along the oriented edge I (its generalized derivative is the cur- 
rent). In [3], the authors have proved the following result (which trivially leads to 
a fluctuation theorem for large deviations), that we state here for a general initial 
distribution (only steady states are considered in [3]): 

Fact 3.3. For any initial distribution v and for any Ai, . . . , A s € R, the limit 

Q(Ai,...,A s ) := lim — lnE„ [e~^=i (3.8) 
t— >oo t L 

exists, is finite and does not depend on v. Moreover, it holds 

Q(Ai, . . . , A s ) = Q(A(d) - Ai, . . . , A{C S ) - \ s ) . (3.9) 

The above Fact 13.31 restricted to the steady state corresponds to formula (39) 
in [3]. As observed there, under Schnakenberg's conditions (see the Introduction), 
this formula leads immediately to the fluctuation theorem for currents restated in 
terms of macroscopic forces (see page 124 in [3]). In Section [6] we show that Fact 
13.31 is a simple corollary of Fact 13.11 This was already pointed out in [3] , by different 
arguments. 
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3.2. Fluctuation theorem for generalized currents associated to a basis of 
oriented cycles. In this section, we extend the fluctuation theorem to currents 
referred to (what we call) a basis of oriented cycles in Q. This answers a natural 
conceptual question. In addition, in applications, one can have identified a nice basis 
of oriented cycles and desire to work with it, without looking for a nice fundamental 
set of oriented cycles (see Section H]). First we fix some general concepts, applying 
some ideas of linear algebra to oriented cycles. 

We say that the family of oriented cycles C\ , . . . , Ck in Q generates all the oriented 
cycles if for each oriented cycle C there exist real numbers a±, . . . , at satisfying C = 
X)i=i a iCi (recall Definition l3.2j) . Such a generating set is called basis if it is minimal, 
in the sense that it does not contain a smaller generating subfamily. It is simple to 
check that minimality can be replaced by independence, i.e. there are not constants 
ai, . . . ,0^ (not all zero) such that Yli=i a £i = where denotes the degenerate 
cycle with S^) = for each oriented edge t. In the last case, we also say that 
C\, . . . , Cfe are independent. Lemma 17.11 in the integrating Section [7] will report 
some very intuitive facts concerning bases, generating sets and independent sets. 
Here we recall that all bases have cardinality given by the chord number s, that 
any fundamental set of oriented cycles is a basis (the opposite implication is false, 
see Section [7J and that for each oriented cycle C the coefficients a\ , . . . , a s in the 
decomposition 



where C±, . . . , C s is a basis, are univocally determined. If the basis is not a fundamen- 
tal set of oriented cycles the coefficients a\ , . . . , a s do not have a simple geometric 
characterization as in (|3 . 5[) . As showed in Section if C\ , . . . , C s is a basis and 
C[, . . . ,C' S is a fundamental set of oriented cycles with associated chords £±,...,£ s , 
then it holds 



(ai, . . . ,a s ) = (Se 1 (C), . . . , S£ S (C))B 1 , B = {Bij)i<ij< s , Bij = 5^.(Cj) . 



Let us come back to our random walk on the graph Q. We want to associate to each 
trajectory up to time t an oriented cycle Ct- We do it as follows. Of the trajectory 
up to time t we record the sequence of oriented bonds (&i, . . . , b n ), hi G £ Q , which 
the walker moves along, one after the other (n is the number of jumps performed 
up to time t). We denote this path as reduced trajectory. Recall the definition of 
the function Ni(t) given after (|3.2p . Note that i\fy(t) depends only on the reduced 
trajectory, since iV^(t) = J2?=i[K^i = ~ ^(^ = Given two different states c, a' 
we fix a path 7 CT)0 -' in Q from a to a' specifying only the visited states and the edges 
along which the path evolves (there is no knowledge of jump times), i.e. 7o- )(T ' is 
represented by a sequence of oriented bonds. We then introduce a random oriented 
cycle Ct on Q as follows: 

Definition 3.4. Fixed a family {'y (Tt<J i} a -/ :< ji , let (bi,...,b n ) be the reduced trajec- 
tory from Xq to X t . If Xq = X t we define Ct as the oriented cycle given by the 



.s 




(3.10) 



(3.11) 
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reduced trajectory itself, otherwise we define Ct as the cycle (b±, . . . , b n , ci, . . . , cv) 
where 7x ,x t = (ci, . . . , c r ). 

We point out that 

N e (t) = S t (Ct) + O(l) W e £ (3.12) 

where the errors 0(1) are uniformly bounded as £ varies in £ Q and t varies in [0, oo). 

From now on, we refer to a basis Ci, . . . ,C S and a family {'y aa i} CTf L a i , fixed once 
and for all. Due to the above discussion, we know that the random numbers 
ai(t), . . . ,a s (t) such that 

s 

Ct = ^ai(t)Ci, (3.13) 

i=l 

exist and are univocally determined. Note that these random numbers depend not 
only from the basis but also on the paths j a y- On the other hand, choosing other 
paths ja-y would change the random numbers ai(t) of quantities uniformly bounded 
as t varies in [0, oo) (as a byproduct of (|3.1ip and (|3.12p ). We point out that, by 
(|3.5p and f)3. 12j) , if our basis is a fundamental set of oriented cycles with associated 
chords ii, . . . ,£ s , then it holds ai{t) = A^. (t) + 0(l). In particular, a part an error of 
order O(l), the number ai(t) is simply the flux along the chord li, which equals the 
time-integrated current along £j. Due to this special case, we call ai(t) generalized 
time-integrated current. 

We can finally state our fluctuation theorem: 

Theorem 3.5. Fix a basis C\, . . . ,C S and paths {j a ,a'} cr ^ (7 i ■ Then, for any initial 
distribution v and for any Ai, . . . , A s £l, the limit 

Q 6 (Ai,...,A s ) := lim -ilnE„ [e-^^l (3.14) 
exists, is finite and does not depend on v. Moreover, it holds 

Q 6 (Ai, . . . , X s ) = Q b {A{Ci) - Ai, . . . ,A(C S ) - X s ) . (3.15) 

The index b refers to the term "basis". Due to the fact that ai{t) = N^.{t) + 0(1) 
when the basis is a fundamental set of oriented cycles, the above result extends Fact 
13.31 The proof of Theorem 13.51 is given in Section 17.11 



4. Applications to molecular motors moving along a polymer 

In this section we consider a molecular motor p3] moving along a polymer (math- 
ematically, a one dimensional periodic environment), under the effect of an external 
force /, using chemical energy to produce mechanical work. We write Afj, for the 
chemical potential difference associated to ATP hydrolysis. We denote by xt the 
position of the molecular motor along the filament and by zt the number of hy- 
drolyzed ATPs minus the number of synthesized ATPs such that zq = 0. We want 
to apply Theorem 13.51 to prove in full generality the relation (jl.2p described in the 
Introduction. 
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The above fluctuation relation (|l,2p has been proved in [2H [22] for a special model 
with two chemical states. The method followed there relies on the manipulation 
of 2 x 2 matrices (since only two chemical states are considered there). On the 
other hand, keeping in mind the above observations concerning macroscopic forces 
and affinities, one would expect the above symmetry (|1.2p to be universal. This is 
indeed what we prove below as confirmation of the large flexibility of the fluctuation 
theorem referred to oriented cycles. To deal this problem in full generality one has 
to work with a very large graph. 

First of all we need to fix a discrete kinetic model. The natural modeling of 
molecular motors is by continuous models called ratchet models [161 128] . Methods 
to derive a discrete model from the continuous ratchet model have been developed 
in full generality (cf. [31] and reference therein), thus leading to the following class 
of kinetic models. The model is a continuous-time Markov chain with state space 
ZxT,r being the finite set of the chemical states of the motor, with master equation 

d t p t {x, a) = uj a x _ hx p t {x - 1, a) + uj a x+l ^ x p t (x + 1, a) - [u%^_ x + uj^ x+1 ]p t {x, a) 

Above, pt(x,a) denotes the probability of the motor to be at site x5 (5 being the 
length of spatial unity in the discretization) at time t in chemical state a. u) xx i is 
the probability rate for a jump from the position x5 to x'5 (with \x — x'\ = 1) if the 
chemical state is a , and , is the probability rate for a chemical transition from 
a to a' if the mechanical coordinate is x. When considering motor proteins moving 
along polymeric filaments, the above rates are periodic function in x, with the same 
period. 

In order to keep information about the ATP consumption, it is convenient to en- 
rich the above discrete model, extending the state space and distinguishing between 
active and passive chemical transitions (i.e. related to hydrolysis/synthesis of ATP 
or to thermal noise) [27]. The new state of the system (motor plus environment) 
is now described by the triple (x,a,z) eZxTx Z, where z denotes the algebraic 
number of hydrolyzed ATP molecules (in the sense that z increases of one unity for 
each ATP hydrolysis and decreases of one unity for each ATP synthesis). We write 
w^'L for the probability rate of a chemical transition from a to a' at the mechanical 

state x with the consumption of I = —1,0,1 ATP molecules. Hence, uf2 , is the 
probability rate for the jump (x, a, z) — > (x, a' , z + I). The other possible jumps are 
the mechanical ones (x, a, z) — > (x r , a, z), \x — x'\ = 1, having probability rate oj x x ,. 
The resulting master equation is the following: 

d t p t {x, a, z) = u)l_ l x p t (x-l, a, z)+u^ +1>x p t (x+l, a, z)- [u) (T XtX _ 1 +u: a x , x+l ]p t {x, a, z) 

+ i u ^,aPt( x ' CJ '' z-l)- uj*' l a ,p t (x, a, z)] . (4.2) 

£=-1,0,1 a'^a 

In what follows, we suppose that all the rates are positive. In principle, this is not a 
restriction since positive but very small rates correspond to very unlikely transitions. 
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Recall that / denotes the load force and A/i the chemical potential difference 
associated to ATP hydrolysis. We write V a for the potential energy of the molecular 
motor in the chemical state a due to the interaction with the polymer. Again, 
the functions V a , a G T, are periodic functions of the same period. The energy 
associated to the state (x, a, z) is given by 

E(x, a, z) = V a (x5) — fx5 — zA/i + const. (4-3) 

Then, the detailed balance condition reads 



x,x+l 



/u% +ljX = exp {-p[V a (x6 + 5)- V a (xS)] + (38 f} 



= exp {-P[V a ,(x5) - V a (xS)] + piAfj,} 
We can now formulate our result: 



(4.4) 



Fact 4.1. Recall that xt, Zt denote the 'position at time t of the molecular motor 
and the total ATP consumption up to time t. Denote by f(3 and xt the adimensional 
quantities ff3 x [lm] and Xt/[lm] ([lm] being the length unit). Then for any given 
initial configuration (xo,Zo = 0,o"o) and for any constants A, 7 G R the following 
limit exists: 

0(A )7 ) :=/nn -ilnE (:COi2o=0i(To) [ e - A ^-^] . (4.5) 

Moreover, it holds 

tf(A, 7 )=tf(7£-A,/3A M - 7 ), A,7 G R. (4.6) 

We point out that the model presented in [T5J, [2TJ E2] is indeed a coarse-graining 
approximation of the model we have described above when considering two chemical 
states and suitable potentials. 

4.1. Proof of Fact 14.11 In order to apply Theorem 13.51 we need to work with a 
finite state space. This is obtained by a suitable projection. 

Recall that the filament is periodic. We define N as the number of 5-units con- 
tained in a spatial period. This implies in particular 

= VJxS + N5), u*' 1 , = cj x +/' 1 , , x 

J x+l,x — ^x+N+l^+N ' u x,x+l = ^x+N^+N+l ■ 

We write for the quotient space TLjNTL given by integers modulo N. Of course, 
sums in Z^r are modulo N. Due to (I4.7D both the potentials and the jump rates can 
be thought of with the spacial parameter x varying in Z^v- 

We introduce a finite graph Q with multiple edges and no loop, indicating the 
canonical orientation of edges. To this aim it is convenient to label the states in T 
as 01, (T2, • • • , <7 m (m = |r|). Then, the vertexes of Q are given by the pairs (x, <7j), 
x G Zjv and 1 < i < m. We put an oriented edge e^* x+1 from state (x, o~i) to state 
(x + 1, <7j). Moreover, given 1 < i < j < m, we put three oriented edges from state 
(x,(7j) to state (x,aj) labeled by the index I = —1,0,1 and we call them e%^ aj . 
See figure [TJ Note that Q has Nm vertexes and Nm + 3Nm(m — l)/2 edges. In 
particular, its chord number is 

s = 3Nm(m- l)/2 + 1. (4.8) 
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(i,(Ti) ~(x + l,<Ti) 

Figure 1. The edges (Q x+V e%£, <£? ><TJ , e% 



_x,l 



We now consider the continuous time random walk Xt on the graph Q, where the 
probability rate for a jump along the edge e u xx+l equals uj x * x+1 , while the probability- 
rate of the reversed jump is oj^+i x - Similarly, given 1 < i < j < m and I = 

— 1, 0, 1, the probability rate for a jump along the edge e%^ aj equals w^,^ , while the 

probability rate of the reversed jump is uj^ T '~l i . Recall that, given an oriented edge 
£, we write Ng(t) for the number of times the random walk jumps along the edge I 
minus the number of times the random walk jumps along the reversed edge £, up to 
time t. Trivially, the above random walk Xt can be coupled with the Markov chain 
considered in the previous subsection getting the identities 



x t /8 



E E ^(0 

x&n l<i<m 



E E K 



(t)-N,.- } (t) 



This fact, together with (|3.3p . implies that the limit (|4.5p exists and is finite, hence 
the function #(A, 7) is well defined. 

It remains to prove (|4.6p . To this aim recall how we have associated the cycle 
Ct to the trajectory of the random walk up to time t (see Definition 13. 4p . Trivially, 
a part an error of order one uniformly in t € [0, 00), the above identities can be 
rewritten as 



Xt/S = Y.xez N Ei<i< m S e °i , , (°t) + 0(1) 



Y.x&L N El<i<7<m \ S t£,± B . ( C <) - S e%'-\ + ■ 



(4.9) 



We give now two different applications of Theorem [33J In the first case we exhibit a 
nice basis which is not a fundamental set, in the second case we exhibit a spanning 
tree leading to a nice fundamental set. 

4.1.1. A nice basis. We introduce a basis of oriented cycles on Q. We call Cq the 
cycle described by the ordered family of oriented edges 



Given x € 

ED 



^0 — ( e l,2i e 2,3> • • • ' e N-l,N' e N,l) ■ 

jjv and 1 < i < j < m, we introduce the three oriented cycles (see figure 



1 x,0 p x ~ l \ 



1X,1 



( x,l -x,0 \ 
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<x,x+l _ I oi z+1,0 ffi -x,0 \ 

fi,aj ri.i+l'^i.ffj ' c a;,a;+l' c cr i ,o-j/ • 




(x,ai) (x + l,(Ji) 

Figure 2. The cycles C^~), Cp* aj , Cp^ 1 . 



Consider now, for each x, the complete graph on the set {(x,o~i) : 1 < i < m} 
with canonically oriented edges e^ iCTj , 1 < i < j < m (see figure [3]). Trivially, the 
ed ges &cr\,cF2 1 60-2,0-3 , • • • ,Com„i,o m form a spanning tree. We call Cal,a-j the oriented 
cycle associated to the chord e%f^j , with (i, j) {(1, 2), (2, 3), . . . , (m — 1, m)}. The 
orientation of Cp^ ai agrees with the one of the chord eZf, aj ■ 

0,04) (x, £T 3 ) 




(x,cri) (^,o- 2 ) 
Figure 3. The complete graph with edges eZf^j, 1 < i < j < m := 4. 

Lemma 4.2. Consider the family of oriented cycles 

(i) Co,_ ] 

(ii) Cct-.oj urai/i x £ Zat and I < i < j < m, 

(iii) Cct-^oj u«f/i x € Ztv 1 < z < j ' < 

(iv) Cct-°oj with x e Zjy, 1 < i < j < m and (i,j) g {(1,2), (2,3), . . . , (m — l,m)}, 

(v) Caf t a i+1 with x € Zjv 1 < i < m — 1. 
ITien t/ie above family of oriented cycles is a basis. 

Proof. The oriented cycles of type (i),(ii),...,(v) are respectively 1, Nm(m — l)/2, 
Nm(m — l)/2, Nm(m — l)/2 — iV(m — 1), iV(m — 1). It follows then that the above 
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family of oriented cycles has cardinality 2>Nm{m — l)/2 + 1, which is the chord 
number s given in ([4. 8ft . 

By point (v) in Lemma [7.14 in order to show that it is a basis it is enough to prove 
independence. Suppose to have a linear combination of the above oriented cycles 
which equals the degenerate cycle 0. First we observe that the edge e£-^ belongs 
only to Cal~a), the reversed edge does not belong to any oriented cycle in the family. 
By applying the operator St with I = e^^*, and then varying x,i,j, we get that 
the cycles of type (ii) do not appear in the linear combination. The same holds for 
cycles of type (hi) (think to the edges e%' ii(r A. Keeping in mind this consideration, the 
edge e x a f >aj with x 6 Zn, 1 < i < j < m and {(1,2), (2,3), ... ,(m- l,m)}, 
appears only in the oriented cycle CpP^ of type (iv) while its reversed edge does 
not appear in any cycle of the family (recall the properties of a fundamental set of 
oriented cycles). Note that edges of the forms e%\ y(Tj or e%\ )(Tj appear in the oriented 
cycles of type (v) but it must be j = i + 1. 

At this point, we know that the linear combination involves only Co and the 
oriented cycles of type (v), i.e. we have 

«oC + £ E 0- ( 4 .io) 

x£Zff l<i<m 

Applying to both members the operator Si with t = e°* x , 1 (i € Zjv and 1 < i < m), 
we get for any x € 7L^ 

alfX+i = o^mii+a ifl<i<m-2, (4.11) 

X,X + 1 _ 

The above system trivially implies that all coefhcients in (|4,10p are zero. Hence, the 
oriented cycles in the above family are independent. □ 

We point out that the above basis is not a fundamental set of oriented cycles. 
Indeed, each edge of Co belongs also to some oriented cycle of type (v), hence Co can 
not contain any chord. 

Let us now compute the affinities associated to the oriented cycles in our basis. 
Due to the detailed balance relations ()4.4p . we have the following weights: 



w{e% iX+1 ) = -P[V a (x5 + 5)- V a (xS)) + (35 f , 



(4.12) 



Hence, 



A{C ) = PNf5 , .4(C'; : .J ) = /3A/i , A{C*£ aj ) = /3A/x , (4.13) 
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(note that conservative force fields never appear in the cycle affinities). We write 



C t = a (t)C + E E 



,x,-l , 



E E 



E E Cw c ^r ( 4 - 14 ) 



\<i<j<m 
(i,j)£{(l,2),..,(m-l,m)} 



xGZjv l<i<m 



By applying (|4.9p to the above decomposition and setting 5 = S/[lm], we get that 
Xx t + jz t = XNSa (t)+ 1 E E ® + <%U (*)] + ^ • ( 415 ) 

The above identity together with Theorem 13.51 implies that 



0(A, 7 ) = Wo -> AA5 , C£~ 7 , C*£ 



7 



x,x+l 



0) 



Q fe (C -> PNf5-\N6 , C^-J -> , C*-^. -> /3A/x- 7 , C*&.. , (^J 1 -»> 0) 

W-7,/3A/i- 7 ). (4.16) 

This concludes the proof of Fact 14.61 



4.1.2. A nice fundamental set. The above arguments can be applied also to funda- 
mental sets of oriented cycles. As example of spanning tree, one can take the tree 
given by the edges with 2 < i < m and x € Z^v, and e a x 1 x+1 with 1 < x < N — 1 
(see figure Alternatively, one can take as spanning tree the one given by the edges 
e CT i,<jj with € {(1, 2), (2, 3), . . . ,(m - l,m)}, and e^ x+1 with 1 < i < m and 

1 < x < N — 1. In both cases, the associated fundamental set is tractable. The 
advantage here is that one does not have to check that the associated fundamental 
set is a basis, since this is automatically. We leave the details to the interested 
reader. 



(1,^4) 



V\ (1,0-2) 































X 



1 



Figure 4. Example of spanning tree where N = 3, m = 4. The 
spanning tree is in boldface, edges e-a'^aj have been omitted. 
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5. Mathematical integration to Section [2] 

The decomposition of the generator L in terms of the rates k^(a, a') described in 
Section [2] is naturally associated to the following construction of the Markov chain 
Xt. First we observe that the evolution of the Markov chain is univocally determined 
by the sequence o"o, o\ , oi , . . . of the visited states (listed in chronological order) and 
by the times T n , n € N, of the n-th jump (with the convention To = 0). In order 
to specify these last quantities, we consider a family of independent exponential 
random variables T^'^\ parameterized by the integers n € N and i : 1 < i < m 
and by the ordered pairs (a, a') of distinct elements in S. The random variables 
T^fr are all defined on a same probability space (fi, J 7 , P), T^fP is an exponential 

variable with parameter (i.e. 1/mean) k^(a,a' ). Note that if'f = oo if kf, = 0. 

The Markov chain Xt starting in <to can be inductively constructed as follows as 
function on the space 0: known the n-th jump time T n and the n-th visited state 
a n , we set = m.m i>a / T^ n J,. Then T n+ \ := T n + T^ n \ i.e. after the n-th jump 
the Markov chain remains at a n for a time , after that it jumps to the state a' 
such that for some i it holds T (n) = T} i,n J,. As well known (cf. [26] [Theorem 2.3.3]), 
the index i and the state a' such that = ']~( l ' n \ are univocally determined a.s. 

Recall the interpretation of the Markov chain as random walk on the graph Q, 
given in Section [2j Known a n and T n , we can think that the walker moves at time 
T n+ \ = T n + T^ to the state a' along the z-edge if = T^'^,. In addition, recall 
that in Section [21 given the initial distribution v on S, we have defined as the 
probability measure describing the evolution of the Markov chain keeping knowledge 
of the crossed edges. We point out that P u is a probability measure on the space fl, 
where all the exponential variables T^'?^ are defined. 



5.1. Proof of Theorem 12.11 The proof follows the strategy stated in [24] [Section 
2.3], i.e. one has simply to extend the proof given in [24] [Section 2.1] for m = 1. 
The only difference is that one has to keep in mind the special construction of the 
Markov chain Xt given above. We give the proof for completeness and also to show 
that the assumptions required in [24] [Section 2.3] are redundant. 



We fix Ai, . . . , A m and set gt(cr) = K a 



It is well known that the 



function [0, oo) 3 t — > gt(cr) £ (0, oo) is differentiable. We think of Xt as a random 
walk on the graph Q with multiple edges and no loop as explained in Section [2j 
Fixed e > we introduce the following events: A is the event that in the time 
interval [0, e] the walker does not jump, B®, is the event that in the time interval 
[0, e] the walker makes a unique jump and this jump takes place along the i-labelled 

(i) 

edge from a to a , while C denotes the event that neither A nor any event B a a , 
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takes place. By the Markov property we can write 



m 

+ y, f2 E °[ e ~ EJLlXjWU){£) 



a':a'^a 1=1 



(i) 



+ E CT 



g t (a)F a (A)+ £ £ e - A ^T CT (i?W^ 5t (a')+E CT 

o':a'±a i=l 



Using that P CT ( A) = l-r(a)e+o(e), IPo-(-B^i') = ^) e +°( £ ) and = o( e )» 

where 

r(a):= £ *(*,</) = £ £ °0 , 

a'-.a'^a a'.a'^a i=l 

the above expansion (|5.ip implies that 

^( fJ )= E [ E ^)(a, ( 7') 1 ^fc (4) (^^) A! ]5t(^)-K^(^)=:(^)(^)- 

a'-.a'^a i:l<i<m 
fc( i '(o-,o-')>0 

The above differential equation and the fact that go(cr) = 1 for any a imply that 
g t (a) = (e tc ±go)(cr) = Y,a>es[ etC -}aa'- In particular, we can write 



e-^MW^) gt{ae) , C 



(5.1) 



E, 



(5.2) 



where (•,•) denotes the Euclidean scalar product in M 5 and 1 denotes the vector 
with all entries equal to 1. 

We point out that the C\ has off-diagonal nonnegative entries. In particular, for 
a large enough, M := C\ + al has nonnegative entries and one can apply to it the 
Perron-Frobenious theorem (cf. |15| [Theorem 1.5.4]). From the fact that Xt is an 
irreducible Markov chain, we derive that C\ is irreducible (i.e. given a ^ a' there 
exists a path a = <ro, <ti, . . . , a n = a' such that [Cx\ aj>aj+1 > for all j : < j < n) 
and as a consequence the same conclusion holds for M. Applying now the Perron- 
Frobenius theorem we conclude that M has a maximal eigenvalue which is simple 
and is associated to an eigenvector v having all positive entries. Hence, the same 
conclusion holds for C\. We call — e(A) the maximal eigenvalue of C\ associated to 
v. We write Vj for the j-th entry of v and set v max = max-,- vj, u m ; n = min-,- Vj. Since 
both the vectors v, v and the matrix e ^ = e~ ta e tM has only nonnegative entries, 
we conclude that 

v-l x e-^(u,v)=v-l x (u,e^v) < (u,e^l) = E„ [ e -E£xW«W 

<v-l n (v,e^v)=v^ n e-^(v,v). (5.3) 
Since {u, v) > the above estimate trivially implies (|2.5[) . 
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Since £(i_ Al i-\ m ) = £(Ai,...,A ro ) (the l.h.s. denotes the transposed matrix, recall 
the definition of C\) and since 1-A \ and £(i-Ai,...,i-A m ) have the same 

eigenvalues, we immediately get (|2.6|) . □ 

6. Derivation of Fact 13.31 from Theorem 12.11 

Recall that we have labeled the canonical oriented edges in £ c as £i,£2, ■ ■ ■ ,£ s , 
£ s +i, ■ ■ ■ An, where £\, . . . , £ s are the chords. Due to Definition ^. 41 the representation 
PTE?]) and due to (EHJ) we get 



q({\ e : I G £ c }) := Urn lnE„ 



lim In E t 

t— ¥00 t 



.1) 



Due to (|3.5p (recall the definition (|3.2p of the weight w{£)) we can write 

s s n n 

£A i s A (c t ) = 2>s A (c t ) + E "(MSttiCt) - E ™&)Sk{Ct) 

i=l i=l i=s+l i=s+l 

s n n s 

= ^Ai%(c t )+ ^)SiACt)- E Y, w ^ s ^ s ^ 

i=l i=s+l i=s+l j=l 

s n n 

= E^" E <ta)S ia (C i )]S ii (Ct)+ E ™(^(Ct). 

i=l a=s+l i=s+l 

The above identity and (|6.ip imply that the limit (|3.8p exists, is finite and equals 



Q(Ai, . . . , A s ) = q({\i - E w(£ a )Si a (Ci)}i<i< s ,{w(£i)} s+1 <i< n 

a=s+l 

Due to the fluctuation relation (|3.4p . the expression in the r.h.s. equals 

n 

l<i<n ) • (6.2) 

a=s+l 

As already observed, given i : 1 < i < s, SV (Cj) = 5ij for 1 < i,j < s. In 
particular, w(£i) + ^2a=s+i w (^a)Se a (Ci) simply equals the affinity A{Ci). Combining 
this observation with the above identities, we get that 



HOD = q[{A(d) - Xi}i<i< s ,{0} s+ i<i<n) = Q({A(d) - Xi}i<i< s ) . 
This completes the proof of Fact 13.31 □ 
7. Mathematical integration to Subsection 13.21 
We list some simple properties concerning linear combinations of oriented cycles: 

Lemma 7.1. The following holds: 

(i) any generating set of oriented cycles has cardinality at least the chord number s, 

(ii) any independent set of oriented cycles has cardinality at most s, 
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(iii) any basis has cardinality s, 

(iv) any generating set of oriented cycles of cardinality s is a basis, 

(v) any independent set of oriented cycles of cardinality s is a basis, 

(vi) given a basis C±, . . . , C s and given an oriented cycle C, the coefficients a%, . . . ,a, 
such that C = X^i=i a i^i are univocally determined, 

(vii) any fundamental set of oriented cycles is a basis. 

Proof. We first prove (i), taking a generating set of oriented cycles C\, . . . ,C^. Then 
we fix a fundamental set of oriented cycles C[ , . . . , C' s and call i\ , . . . , £ s the associated 
chords. We introduce the k x s matrix B with Bij = (Ci). Due to (|3.5f) it holds 
Ci = ^2j = i BijC'j for all i : 1 < i < k. On the other hand, by definition of generating 

set, there exists a s x k matrix A such that Cj = 5Zi=i ^-ji^i for all j : 1 < j < s. 
Then, for all r : 1 < r < s, it holds 

k k 
8j tr = Si r (Cj) = ^2 AjiSi r (Ci) = ^2 AjiBi r . 

i=l i=l 

The above relations can be rewritten as I = AB, where I denotes the s x s identity 
matrix. This trivially implies that B must have rank at least k, hence k > s. 

In order to prove (ii) suppose that C±, . . . is an independent set of oriented 
cycles. Taking a fundamental set C[,...,C' S as above, for a suitable k x s matrix 
we can write Ci = J2j=i BijCj for all i : 1 < i < k. If k > s, then we could find 
a non trivial zero linear combination of the rows of B, i.e. not all zero coefficients 
a\,...,ak such that YM=i a i^ij = for all j : 1 < j < s. Since Ya=i o-iSi{Ci) = 
Yli=i ^j=i a iBi,jSi{C'j), this fact trivially implies that ^i=i a «^« = m contradic- 
tion with the hypothesis of independence. This completes the proof of point (ii). 
We can directly prove point (v). To this aim suppose that k = s. Then, as ar- 
gued above, independence implies that the above matrix B has rank s, hence it is 
invertible. Calling A := B , since Ci = Y2j=i ^ijC'j , this trivially implies that 
X)f=i A-kiCi = C' k , for all k : 1 < k < s. In particular, the family C\, ■ ■ ■ ,Ck generates 
a fundamental set of oriented cycles which generates all oriented cycles by (I3.5p . 

Point (iii) is an immediate consequence of (i) and (ii). To prove point (iv), let 
C\, ... C s be a generating set of cardinality s. If these cycles were dependent, then 
we could find a generating subset of oriented cycles, in contradiction with point (i). 
Property (vi) follows from the definition of basis, while property (vii) follows from 
(iv) and d£5J). □ 

We point out that not any basis is a fundamental set of oriented cycles. Consider 
figure 5. Vertexes are numbered from 1 to 6. Consider the oriented cycles C\ = 
l->-2-»5->-6-»l, C 2 = 2-»-3-»4-»-5-»2, C 3 = 3-»4->-6->l->-3 
and C4 = 4— >5— s>6— s>4 (above we have indicated the visited vertexes, in order of 
visit). This family is a basis. Indeed, it has cardinality given by the chord number 
and moreover it trivially generates the fundamental set of oriented cycles associated 
to the maximal tree given by the edges {1,2}, {1,3}, {1,6}, {4,6}, {5,6} depicted 
in the figure by boldface (with exception of the edge {1,3}). On the other hand 
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Ci , . . . , C4 do not form a fundamental set of oriented cycles since each edge of C4 
belongs (with the same or with opposite orientation) to some other oriented cycle 
Ci, i 7^ 4. Hence, no edge of C4 could be the chord associated to a hypothetical 
maximal tree. 




2 



Figure 5. Graph for an example of basis which is not a fundamental set. 



We conclude this section with the proof of Theorem [37 



7.1. Proof of Theorem 13. 5L We fix a fundamental set of oriented cycles C[, . . . ,C' S 
and call l\, . . . ,t s the associated chords. Recall that the orientations of C[ and li 
agree. We introduce the s x s matrix B defined as Bij = 5^.(Cj). Then by (|3.5[) 
it holds Cj = X^i=i BjiC'i- We point out that the matrix B is invertible (see the 



arguments used in the proof of Lemma 17.11) . 

Given an oriented cycle C we can express the coefficients in the representation 
C = J2t=i a iCi by means of the fundamental set of oriented cycles and associated 
chords as follows. Applying the operator S^. to the identity C = Y2t=i a iCi-> we S e t 

s s 

S £j( C ) = '^2 a iSe j (Ci) = ^2aiB {j , j = l,...,s. (7.1) 

i=l i=l 

This implies that 

(a 1 ,...,a s ) = {S ll {C),...,S ls {C))B- 1 . (7.2) 

Writing < •, • > for the scalar product in M s and denoting by S(t) the column 
vector with entries S^ifit), ■ ■ ■ , Se s (Ct), the above equation (j7.2|) applied to the ori- 
ented cycle Ct gives ^i=i a «(0^« = < S{t),B~ l \ >. In particular, the expectation 
in (|3.14p can be written as 



-<S(t),B- 1 A> 



By the above identity, the representation (I3.12P and Fact I3.3[ we get that the limit 
in (I3.14p exists and it holds Qb(^) = Q(-B -1 A). Applying the fluctuation relation 
(|3.9p and afterwards again the identity Q b {-) = QiB^ 1 -), we get 

Q6(A) = Q(B _1 A) = Q(A - B^X) = Q b {B(A - B _1 A)) = Q b (BA - A) , 
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where A denotes the column vector with entries the affinities A(C[), . . . , A(C' S ). To 
conclude, it is enough to note that 



3=1 



s s s 

(BA)i = £ ByAify = Si, (QAiC'j) = £ S tj (CO [ S e (fyw(£) = 

j=i 3=1 ie£c 

s 

£ w(£)[ St, (QSeity] = Y,w{l)S i {C i )=A{C l ). 



3=1 
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